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I. INTRODUCTION 

In circuit QED the model of a two-level system in- 
teracting with one or more harmonic oscillators can be 
implemented combining a superconducting qubit with a 
microwave resonator or a transmission lineal Compared 
to exp eriments in Quantum Optics with microwave cavi- 
tieiPEl or with trapped ions^, the superconducting circuit 
experiments have one important advantage: the strength 
of the qubit-photon coupling. The fact that supercon- 
ducting resonators and superconducting qubits follow es- 
sentially the same physical laws makes it possible not 
only to reach the strong coupling regime^', in which mul- 
tiple Rabi oscillations are possible within the decoherence 
of the cavity or the qubit, g 3> k, 7, but also entering the 
ultrastrong coupling regime, g ~ w, in whi ch t he inter- 
nal and interaction energies become similaJ^. In this 
new regime the dynamics is very fast and the usual ap- 
proximations such as the Rotating Wave Approximation 
(RWA) in the Jaynes-Cummings model break dowrP^^. 

One of the most astounding predictions of the ultra- 
strong coupling regime is that a single qubit can distort 
its electromagnetic environment, giving rise to a ground 
state in which the qubit is dressed with photons. As we 
will show in the following sections, in the case of a qubit 
and a single harmonic oscillator, this translates into a 
state which is a superposition of a deexcited qubit and a 
vacuum, with other states in which the qubit, the oscilla- 
tor or both are pop ulated with excitations and photons, 
respectivel y 10 * 12 * This is a completely non-RWA effect 
which requires large values of the coupling to be observed. 
More precisely, the excitation probability grows approx- 
imately as p e cx (g/uj) 2 and g has to become comparable 
to the energies of a photon, Hu, or of a qubit, Thuq, making 
the interaction dynamics both very strong and very fast. 
From the experimental point of view it would thus seem 
unfeasible to probe a physics that takes place at speeds of 
uj ~ 1 — 10 GHz, while the typical measurement appara- 
tus in circuit-QED have response times which are much 
slower, of about 50 ns. There are four routes to escape 
this problem: making the ultrastrong coupling switch- 
able by design^, dynamically turning it off by external 
drivings^, engineering faster measurement apparatus or 
looking for new ways to extract information out of slow 



measurement devices. 

In this work we take the slow route, showing that is 
possible to extract valuable information from the fast dy- 
namics of the system with current measurement technolo- 
gies. We will study what happens to an ultrastrongly 
coupled qubit-cavity system when the qubit is subject to 
repeated measurements by a detector with a slow repe- 
tition rate that is only capable of performing weak mea- 
surements of the state of the qubit. The main goal is 
to detect the qubit in its excited state starting from the 
ground state of the system. The first measurement has al- 
ready a small probability of success, as commented in the 
previous paragraph. In case of failure the system is pro- 
jected to a non-equilibrium state which rapidly exhibits 
a dynamics with an oscillatory probability of excitation, 
mainly due to non-RWA transitions from the ground 
state of the qubit \g) to the excited one |e). By means 
of performing repeated measurements, we will show that 
the detector is able to probe these usually considered 
as "virtual" excitations of the qubit and the cavity and 
at the same time reveal information of the interaction 
model. More precisely, the repeated measurements ac- 
cumulate information exponentially fast and behave like 
an anti-Zeno effeclP^in which the qubit is projected onto 
its excited state, revealing those ground-state excitations 
that we were looking for. We show that this anti-Zeno 
"decay" \g) — > \e) , is very efficient and does only require 
a short number of repeated measurements with a repe- 
tition rate which is much slower than in the standard 
anti-Zeno effect. 

Like previous proposals for probing the ultrastrong 
coupling limilP^BS the anti-Zeno dynamics in this work 
is supported by the counter-rotating terms in the qubit- 
resonator interaction, using as seed the ground state ex- 
citations of these systems. The phenomenon is absent 
in the limit of RWA in Jaynes-Cummings models. Let 
us remark that the non-RWA effects are being exten- 
sively studied not only in the ultrastrong coupling regime 
of circuit-QED but also in other fields like Quantum 
Optics 19 . Models of repeated measurements on super- 
conducting qubits were considered for insta nce in RefP3 
and have been implemented in the lab 1 21 * 22 !. 

The structure of the text is as follows. In Sect. II we 
will show that the eigenstates of the hamiltonian, and in 



2 



particular the ground state of a qubit-cavity system in 
the ultrastrong coupling regime are not separable, \g, 2n) 
or \e, 2n + 1) , but linear combinations of these vacua and 
excitations. More precisely, the qubit- resonator ground 
state contains a contribution of |e, 1) which grows with 
the coupling strength and becomes relevant in the ultra- 
strong coupling regime, g ~ u. We will see that after a 
few ideal periodic projective measurements of the qubit 
state, the probability of finding that it is in the state 
\g) tends quickly to 0, even if an uncertainty in the time 
taken by the measurement is considered. In section II. 
E we will consider a realistic model of measurement in 
which large amounts of errors are allowed, showing the 
robustness of our method. Section II. F is devoted to the 
analysis of the role of relaxation and dephasing. We con- 
clude in section IV with a summary of our results. 

II. DETECTING SELF-EXCITATIONS OF THE 
VACUUM 

A. The Dicke model 

We will consider the following Hamiltonian, corre- 
sponding to a qubit-cavity system 

H = Hq + gHi = hua^a + ^-a z + hga x {a + a f ), (1) 

where Hloq is the energy splitting between the two levels 
of the qubit |e) and \g) , w the frequency of the photons in 
the cavity or resonator field and g the coupling strength. 

In the weak and strong coupling regimes, in which the 
coupling g <C u), ujq is only compared to the decay rates of 
the cavity and the qubit, one may treat Hi as a small per- 
turbation on top of the bare qubit and resonator states. 
In this limit the counter-rotating terms a*a + : a<j~ aver- 
age out, and the total Hamiltonian becomes equivalent 
to the Jaynes-Cummings model, whose ground state is a 
separable combination of the qubit ground state and a 
cavity vacuum, \g, 0) . 

In this work we are interested however in the ultra- 
strong coupling regime, in which g approaches the qubit 
and photon frequencies, uj and luq. In this case it is more 
convenient to look at the state space in the language of 
parity subspaces^, and treat H and Hj on equal foot- 
ing. Within this picture, the Hilbert space splits up in 
two different chains of states coupled by Hj, and in par- 
ticular the ground state of the system becomes a linear 
combination of states in the even parity sector 

\G) = c \g0) + Cl \el) + c 2 \g2) + c 3 |e3> + . . . (2) 

where the coefficients Cj depend on g, u),u>q. 

B. Detecting excitations with one measurement 

One of the goals of this paper is design a protocol for 
measuring the tiny excitations in the ground states - 
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FIG. 1. (Color online) Probability of excitation for of a qubit 
p e (blue, solid) vs. the dimensionless ratio g/u), for a qubit- 
resonator system [Eq. |TJ] in the ground state of a ultrastrong 
coupled limit 10 = CJo = 1 GHz. This line is undistinguishable 
from a quadratic fit (green, dashed). 

|ci| 2 + |c 3 | 2 + |c 5 | 2 +. . . in Eq. Q — . Let us assume for now 
that we have a good measurement apparatus and that we 
perform a single measurement of the qubit in the ground 
\G) of the system. In Fig. [I] we plot the probability of 
finding the qubit excited after just one measurement 

Pe = (Pe) G =■■ (\e) (e\) G (3) 

against different values of the coupling strength, assum- 
ing always uj = ujq and g/uj < 1. For the strongest cou- 
plings the values of p e are sizable. Moreover, we have: 

2 

Pe=A^, (U=U , ^<1) (4) 

This quadratic behavior comes as no surprise. The main 
contribution to p e is |ci| 2 . If we think of \G) as the 
free vacuum \g, 0) dressed by the interaction Hi, then 
\c\\ may be computed from perturbation theory in in- 
teraction picture, the leading term being proportional to 
|(e, 0)| 2 . It is interesting to see how these con- 

tributions quickly grow as g approaches u, but that at 
the same time the signal in current experiments with 10% 
coupling strengths, might have a too small excitation sig- 
nal to be accurately detected. 

This work is born from the idea that perfect projec- 
tive measurements in c-QED might be too difficult, as 
existing measurement apparatus may be too slow or not 
have enough sensitivity to capture those excitations. The 
constraint of time is found, for instance, in flux qubit 
measurement devices based on SQUIDs, which roughly 
work as follows: A very short current pulse is sent to 
the SQUID, instantaneously changing its potential from 
a periodic function to a washboard potential. In this brief 
period of time, one of the flux qubit states which is sitting 
inside the SQUID may provide, through its intrinsic cur- 
rent and flux, enough additional energy for the SQUID 
to tunnel into a voltage state. This stochastic process is 
random in time and does no have a 100% success rate. 
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Moreover, it requires an additional sustained current that 
keeps the SQUID in that voltage state during an integra- 
tion time large enough for the electronics to realize that 
the measurement succeeded. Adding the excitation and 
integration phases, the best experimental setups bring 
the detection time down to tens of nanoseconds, which 
is still slower than the qubit-resonator dynamics -1.6 ns 
for a 600 MHz coupling, and much shorter for the qubit 
and resonator periods, 

An additional complication of the ultrastrong coupling 
limit is that an arbitrary measurement device might not 
have enough good coupling to either the qubit or the 
resonator in the ultrastrong coupling regime. If we as- 
sume that both quantum systems interact so strongly 
that their eigenstates are highly entangled states with 
large energy gaps, ~ g,(d,u>o> the detector could have 
problems coupling to those states and breaking their en- 
ergy level structure. In other words, the measurement de- 
vice couples through an operator, a z , which typically rep- 
resents a perturbation of the qubit-resonator model, and 
if that perturbation, which aims at breaking the linear 
combinations |2|, is not strong enough, it might not ex- 
tract any information from the system, or the amount of 
information might be reduced, becoming an off-resonant, 
weak dispersive measurement. 

All these considerations brought us to the idea of us- 
ing more than one measurement steps in the same ex- 
periment, with the aim of increasing the amount of in- 
formation that it is extracted from the same state. This 
can be done because the kind of measurements done in 
experiments are non-destructive: the same qubit can be 
continued to be measured at another time. It is true, 
however, that the interval between measurements might 
carry a strong, fast and almost chaotic dynamics [Fig. [2], 
arising from the fact that the measurement brings the 
system into a non-equilibrium state, even if it did not 
produce any information. We will show that this is not 
a limitation, but a plus, and that the repeated measure- 
ments may characterize the intermediate dynamics. 



C. Repeated measurements: survival probability 

If we measure the qubit once, the measurement ap- 
paratus does not click and we are working with a per- 
fect projective measurement, we conclude that the qubit- 
resonator system has been projected onto the state 



I*) = ^2c 2n \g,2n) 



(5) 



which is a (normalized) linear combination of deexcited 
qubits and some photons in the cavity. By measuring 
the ground state in an improper basis, we have created 
a non-stationary state that will evolve very quickly, with 
frequencies that are close to g,u and luq. Lacking any 
other relaxation mechanism than the cavity and qubit 
decoherence times, these oscillations will be sustained for 




FIG. 2. (Color online) After the qubit has been measured 
once, the qubit in the qubit-resonator system is left in a non- 
equilibrium state, | . Here we plot the probability of excita- 
tion for the qubit pi e as a function of the dimensionless time 
ujt, shortly after that measurement. We show three situations, 
ui — cjo = 1 GHz and g/uj = 1/3 (dashed), g/ui — 2/3 (solid) 
and g/to = 1 (dotted), which exhibit fast dynamics. 



a large period of time, causing the qubit to get reexcited 
multiple times. The excitation probability 



Ple(t,0) = (*(«) Pe 



(6) 



may be computed from the initially measured state as 

\*(t))cxe- iHt (l-P e )\G). (7) 

As Fig. [2] shows, p\ e exhibits very fast oscillations, but 
also average to a nonzero value, which is always close 
to the ground state excitation probability of the qubit, 
p e = X)nl c 2n+i| 2 - Consequently, if we perform a sec- 
ond measurement at a later time t\ we will have again a 
certain probability of success pi e (ii,0) of detecting the 
state |e) , and a certain probability of failure 0) = 

1 — pi e (ti,0). In the latter case the system is projected 
to a new state with a new time dependent probability 
P2e{t2,t\), and so on. After a few measurements we can 
define the survival probability as the probability that we 
have never detected a state |e) in the qubit 

Pg = PgPl g (tl,0)p 2 g(t2,tl)...p N g(t N ,t N -l). (8) 

A key idea in the interpretation of this formula is the fact 
that the intermediate probabilities p ng are on average 
very similar, and almost independent of the timespan 
among measurements. For the range of couplings that are 
within intermediate reach in experiments, g/u ~ 0.1 — 1, 
we have verified numerically and perturbatively that this 
probability is well approximated by a quadratic law 



Png 



1 



Xn n 



(9) 



with minor differences among realizations, Xn- The accu- 
mulation of products in Eq. (|S|) leads to an approximately 
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exponential decrease of the survival probability 



N 



exp 



(10) 



as long as << 1- This exponential behavior is typ- 
ical of the so called anti-Zeno effect, in which repeated 
measurements of a quantum system accelerate the tran- 
sition of a quantum system between two states. In our 
case the repeated measurements are rather creating a 
non-unitary evolution that excites the qubit from \g) 
to |e) using as seed the nonzero excitation probability 
Pie = J2 n \ c 2n+i\ 2 which is present in the equilibrium 
state of the qubit-resonator system. This last point is 
particularly important because this anti-Zeno evolution 
is impossible when the ground state of the qubit and 
the resonator is the vacuum \g, 0) . In this case g/co is 
so small, and p\ e so close to zero, that all measurements 
will give no signal at all and the qubit will remain in the 
state | g) for the duration of the experiment. As we will 
see in the following, there is a key difference between the 
effects described in this section and the standard anti- 
Zeno effect: we need only a few measurements and they 
can be widely spaced in time. 

In the following sections we will summarize extensive 
numerical studies of the anti-Zeno dynamics. We have 
contrasted these with various semi-analytical methods, 
one of which, the use of truncated Hilbert spaces, helps 
us in understanding the reason for this behavior. For the 
range of couplings of current interest, g/uj ~ 0.1 — 1, it 
suffices to take two photons, and the ground, |G) , plus 
the two excited states \E'), \E") within the same parity 
subspace. All states can be expanded as in Eq. (pi) with 
coefficients Ci, c^, c"., as linear combinations of \gv), | el) , 
\g2) . After the first measurement, the qubit will end up in 
an excited state with probability |ci| 2 and it will remain 
in the unexcited state with |co| 2 — 1 — |ci| 2 , ending up 
in a combination 



P g \G) = c \G)+c' \E')+ ( %\E* 



(11) 



The crudest approximation would be to neglect all ex- 
cited state contributions and assume that after each mea- 
surement, either the state |e) is detected, or the system 
ends up in \G) . In this case the survival probability would 
be exactly exponential 
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(12) 



In practice, however, the combined system does not end 
up only on the ground state, but gets excited state con- 
tributions from \E') , \E") . When we average the contri- 
butions over the period in which the measurement takes 
place, we find that already after the first measurement 
step, the excited states add up to the total probability, 
(pie) T = \ci\ 2 \c \ 2 + 14 1 2 \c' \ 2 + \ci\ 2 \c'o\ 2 + .. . , enhanc- 
ing the original behavior. 
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FIG. 3. Survival probability after eight measurements Pg vs. 
g/u), with uj = Wo- The measurements are performed with 
periods T\, T 2 = v^Ti and averaged over 100 values of Ti 
within the interval 27r[0.1, 5]. Note how the law approximates 
the Gaussian behavior in Eq. (10 1. 



D. Numerical experiments 

We have verified the anti-Zeno dynamics and the expo- 
nential law ( 10 ) by means of exact numerical simulations 



in which we compute the outcome of repeated measure- 
ments on a qubit-resonator Dicke model ([!]) . We will now 
explain the main results of this study. 

From an experimental point of view it might be inter- 
esting to maximize the exponent x, optimizing the mea- 
surement repetition rate to hit all the maxima in the 
evolution of the excitation probability [See Fig. [2] . How- 
ever we found that this is very difficult and demands a lot 
of precision on the measurement apparatus; for small er- 
rors or some measurement randomization this procedure 
drives the apparatus into exactly the opposite regime: 
always hitting the minima of excitation. Seeking a more 
robust, less demanding approach we opted for using two 
incommensurate periods, Ti and T 2 ~ y/2Ti, simulating 
measurement at times 



t n e {T x , Ti + T 2 , 2Ti + T 2 , 2Ti + 2T 2 , . . .}, 



(13) 



and at most optimizing the value of T\ . 

With this approach, and exploring different values of 
Ti, we have studied the survival probability and con- 
cluded that the exponential laws are really accurate. As 
shown in Fig. [3j if we fix the total number of measure- 
ments to be N — 8 and sample various periods, Ti, we re- 
cover on average the Gaussian behavior exp(— Nxg 2 /to 2 ) 
deduced in Eq. ( 10 ). Instead of fixing the number of mea- 



surements, we can also study the same law and verify the 
exponential decay with respect to N. This is shown in 
Figs. [4^,-b, where we plot the accumulated survival prob- 
ability, P^ , as a function of time, and fit it against the 



same exponential ( 10 ) 



It is important to remark that the exponential decay is 
a robust signature that survives even when the measure- 
ment does not take place at precise times, t n , from the list 
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FIG. 4. (a) Probability after the n— th single measure- 
ment, Pn 9 , (dashed) and accumulated survival probability, 
P™ — ]~[n=i P n 9 ( sou d) vs. dimensionless time cot. We use 
g/ui — 1 and perform measurements with approximate peri- 
ods ojTi = cjT 2 /\/2 = 2tt, averaging over random perturba- 
tions of the actual measurement time, t„, within the inter- 
val ut n + [-0.27T, +0.2tt]. (b) Survival probability (solid 
lines) and the corresponding exponential fits (dashed lines) 
for g/cu = 1/3 (circles), 2/3 (crosses) and 1 (squares), (c) 
Mean value p g (solid) and the corresponding quadratic fit 
(dashed) vs. dimensionless coupling strength All plots as- 
sume to = u>o and (b,c) use wTi = 37r/4 



given before (13). This has been verified by simulating 
multiple runs in which t n is randomly perturbed around 
its average value, and computing the survival probabil- 
ity. We want to remind the reader the importance of this 
robustness, because some measurement apparatus such 
as SQUIDs behave stochastically and produce a signal 
at a random time that can not be determined a-priori. 
The fact that the measurement protocol and the result- 
ing physical behavior are independent of a precise control 
is encouraging. 



The accuracy of the exponential law ( 10 ) suggests that 



the survival probability of a single measurement remains 
constant throughout a single experiment, p g ~ p ng Vn. 
This is qualitatively confirmed by Fig. [4^,, where we show 
that these values oscillate around a mean one that is close 
to the average population of \g) in the ground state, i. e. 
to p g . This suggests us to consider average values and 
approximate 
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(14) 




FIG. 5. Survival probability P* vs. t/T\, with g = uj — ujq 
GHz and ojT\ — ir (blue, circles), 2n (red, crosses) and 
37r(green, squares). Each marker corresponds to a measure- 
ment. 



mate is confirmed by Fig. [4J:, where the quadratic fit is 
almost undistinguishable from the actual behavior. The 
final question which remains to be answered is whether 
the exponent \ depends on the frequency of the mea- 
surements or not. For that we have fixed the coupling 
strength and explored three values of the period, T%, 
studying the average exponential behavior. The result 
is shown in Fig. [5j collapsing all numerical simulations 
in the dimensionless quantity t/T\, and finding that they 
have very similar slopes. 



E. Weak measurements 

So far we have considered ideal projective measure- 
ments, introducing only some stochasticity in the time 
at which the measurement event is produced. We will 
now add another ingredient to our measurement model, 
which is the possibility that the detector only performs a 
partial measurement, leaving the state "untouched" with 
a nonzero probability, e. 

We can easily model an imperfect detector using the 
formalism of completely positive maps, operations that 
transform density matrices into density matrices. If p and 
p' are the states of the qubit-resonator system before and 
after the measurement, we will write, up to normalization 



(l-e)(t-P e )p(±-P e )+ep. 



(15) 



which has the expected quadratic behavior. This esti- 



This is read as follows. With probability e the measure- 
ment device will do nothing, leaving the state untouched. 
With probability (1 — c) the measurement device will de- 
tect the state of the qubit. In this case it will either give 
us a positive signal, moment at which we will stop the 
experiment, or it will not produce anything at all, and we 
will continue with the projected state (1 — P e )p(l — P e ), 
that has the qubit deexcited, \g) . 

This qualita tive model describes measurements from 
a SQUIFJ2il22l, where we place ourselves on the verge of 
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FIG. 6. (Color online) Survival probability (solid lines) 
and p„ g (dashed lines) vs. cut, with g = a) = OJo = 1 GHz , 
uiT\ — 2n and e = 0.2, 0.1 and (crosses, squares and circles, 
respectively). The probability of a measurement at a time 
ujt — ujto is averaged over 20 random values within the interval 
[ut - 0.2-K,uxto + 0.2tt] 



tive terms, an infinite stream of photons leaking out of 
the cavity and other controversial phenomena. 

These effects disappear when one rederives the mas- 
ter equation from first principles, using the qubit-cavity 
eigenstatcs of the ultrastrong coupling model and the 
usual zero temperature baths. In the resulting models 
the main relaxation mechanisms are found to be the de- 
cay to the ground state |G) and a dephasing of the joint 
cavity-qubit states — in other words, dissipation and de- 
coherence in the proper basis — . If we assume this rea- 
sonable model, then we can conclude that the exponen- 
tial laws derived in this manuscript are not significantly 
distorted. To begin with, relaxation to the ground state 
\G) just makes the experiment closer to the truncated 
Hilbert space model considered in Sect. |II C| and in par- 
ticular to the exponential law from Eq. (12 1. For strong 



couplings, decoherence amounts to random modulations 
of the qubit-cavity energy levels, without significantly af- 
fecting the populations, |cj| 2 . Since this is the most rele- 
vant quantity in all the previous discussions, we can also 
expect that, up to minor changes in the rates, the anti- 
Zcno effect will also survive. 



III. CONCLUSIONS 



metastability and assume that if the qubit is in the ex- 
cited state, |e) , the SQUID will tunnel to the voltage 
state with probability (1 — e), giving no signal for \g) . 
Note that with probability e the SQUID may not tunnel 
and then we will gain no information about the qubit or 
the resonator. 

In Fig. [6] we analyze the impact of e in our previous 
results. Even for large errors e = 0.2 we retain the ex- 
ponential behavior observed in Fig. [4ji, with acceptable 
error bars that decrease with increasing number of mea- 
surements — in other words, the qubit is still efficiently 
projected to the excited state. 



F. Relaxation and dephasing 

Throughout this work we have considered in the nu- 
merical simulations the model given by the Hamiltonian 
in Eq. |l]) which do not include effects like relaxation 
or qubit dephasing, usually included in master equation 
approaches. 

We want to remark that it is still an open question, 
both experimentally and theoretically, to understand and 
model the dissipation and decoherence processes of quan- 
tum circuits in the presence of uffrastrong qubit-cavity 
coupling. One popular approac h 1 23 * 24 ! is to combine the 
usual photon leakage mechanism from quantum optics 
models, C(p) ~ 2apa) — a)ap — pa^a, with the qubit- 
cavity Hamiltonian. Note that in such a combination, 
the asymptotic states of the dissipation (the vacuum) 
and of the interaction (populated cavity) are incompati- 
ble, and one may find excitations induced by the dissipa- 



We have considered a system consisting in a supercon- 
ducting qubit coupled to a closed transmission line, op- 
crating in the ultrastrong coupling regime. The ground 
state in such scheme is not just a product of the ground 
states of the qubit and the cavity, as is the case for weaker 
couplings. On the contrary, the vacuum of the system is 
dressed by the interaction and so it contains a relevant 
probability of finding the qubit excited. This probability 
is proportional to the square of the coupling strength. 
We have introduced a protocol for detecting that excita- 
tions with certainty, maximizing the small probabilities 
that are obtained with only one measurement. 

Our main result is that, after a number of periodic 
measurements of the qubit, the probability of finding it 
in the ground state in all the measurements goes expo- 
nentially to zero, even if the measurements are weak and 
are performed with a slow repetition rate in compari- 
son with the fast dynamics of the interaction. We re- 
fer to this as slow quantum anti-Zeno effect. Like the 
well known quantum anti-Zeno effect, the result is the 
acceleration of a transition, in this case the exotic tran- 
sition | g) — > |e) , which becomes relevant in this regime 
due to the breakdown of the RWA. But this procedure is 
less experimentally demanding, since it requires a smaller 
number of measurements and a shorter duration of the 
period at which they are performed. We have shown that 
the protocol is robust to large errors in the measurement 
process, when a realistic SQUID readout is considered. 

This is one of the first experimentally accessible conse- 
quences of the new ultrastrong coupling regime and can 
only be derived beyond the RWA. The physical nature of 
the ground state qubit self-excitations, commonly consid- 
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ered as a virtual process without possible experimental 
record, seems now to be clear. Moreover, although the ul- 
trastrong coupling entails a very fast dynamics, we have 
shown that valuable information of the interaction can be 
extracted efficiently with the current slow and imperfect 
measurement technologies. 

Finally, we want to remark that strong qubit excita- 
tions have also been found theoretically in models that 
combine the full Rabi coupling with traditional dissipa- 
tive contributions^^. However, the form of those dis- 
sipative terms is questionable in non-RWA setups, and 
furthermore, there is no justification to equate the sparse 
measurement setup in this work to a particular dissi- 



pative model. This lack of equivalence between models 
manifests in the fact that, as we have seen numerically, 
the sparsely repeated measurements can hit certain res- 
onances that invalidate the anti-Zeno dynamics. 
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